T radiative decays to light quark jets and color octet mechanism 
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We study radiative decays of T to light quark jets in nonrelativistic QCD by taking both the color 
singlet and color octet bb operators into consideration. The cut for quark jet energy and cut for the 
angle between two quark jets are introduced. The sensitivity to the soft and collinear singularities 
in the loop integrals are greatly reduced by these cuts. With the jet energy cut of about 1 GeV, 
and the jet angle cut of about 36°, the branching ratio for T — * ■yqq is found to be 8.2 x 10 -4 from 
color singlet contributions. The color octet contributions could be much larger than that of color 
singlet, depending on the estimate of the color octet matrix elements. This process may provide a 
new test for the color octet mechanism in nonrelativistic QCD. 
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I. INTRODUCTION 

The observation of T decay into three-gluon jets T — > 
ggg has been viewed as a crucial event for verifying 
quantum chromodynamics (QCD). The radiative decay 
T — * jgg has also been an important testing ground for 
QCD. These processes are computable in perturbative 
QCD due to the large mass of the b quark. Recently 
the CLEO Collaboration has accumulated a large num- 
ber data sample for the low lying T resonances [l|, 
This encourages us to study the radiative decay process 
T — * jqq, where q and q (q — u, d, s) are referred to as 
light quark jets. 

Jets are considered to be footprints of quarks and glu- 
ons produced at short distances. While the gluon jet, its 
production angular distribution, and the determination 
of the gluon spin have recently been studied [3j, less atten- 
tion has been paid to quark jets, say, e.g. in T decays. 
Since the gluon jet is expected to be "fatter" than the 
quark jet[4| or the mean multiplicity of the gluon jet is 
expected to be larger than the quark jet (see e.g. @), it 
will be interesting to measure the quark jets production 
in comparison with the gluon jets in T radiative decays. 

In heavy quarkonium production and decay an impor- 
tant issue is about the color octet mechanism in nonrel- 
ativistic QCD (NRQCD)[6], which is introduced to ex- 
plain the large cross sections of inclusive charmonium 
hadroproduction at large pt in pp collisions measured at 
the Fermilab TevatronQ, and other decay and produc- 
tion processes. Though the color octet mechanism has 
gained a series of successes, some difficulties still remain 
and more tests are needed (for a recent comprehensive 
review, see Q). In this regard, the radiative decay pro- 
cess T — > jqq may provide an interesting test for the 
color-octet mechanism, since the color-octet bb compo- 
nent (with soft gluons) in the T wave function could give 
substantial contributions to this process (for discussions 
on color-octet contributions in T hadronic and radiative 
decays, see, e.g. [1, [l(| Ej], see a ls°@ f° r more refer- 



ences). 

Based on the above considerations, we will study the 
three body decay T — > jqq, where q — u,d, s. We use 
the bb bound state description for T in NRQCD and give 
the results through both color-singlet and color-octet bb 
operators. In NRQCD, with the Fock state expansion 
in terms of v, the relative velocity between b and b, the 
decay width for T(1S) can be written as a sum of contri- 
butions from various bb channels with different color and 
angular momentum: 

dT[T(15)] = dr[66i( 3 5i)]<T|Oi( 3 5i)|T)+dr[66 8 ( 8 5i)] 
x (T|O 8 ( 3 5i)|T) + dr[66 8 ( 1 S , )] 
x (T|O 8 ( 1 5 )|T) + (^(2J+l) 



x dr[6& 8 ( 3 P 7 )])(T|O 8 ( 3 P )|T), 



(1) 



where the long distance matrix elements with different 
color (denoted by the subscript 1(8) for color-singlet 
(octet)) are involved. In Section 2 we will describe our 
method and give a calculation for the color singlet con- 
tribution. In section 3, we will study color octet contri- 
butions. Then we will give a summary in Section 4. 



II. 



COLOR SINGLET CONTRIBUTIONS 



The q and q dijet production through T radiative decay 
at leading order is described by the process 



T(P) - 19*9* -> TfoMWfcs) 



(2) 



as shown in Fig.l (another three diagrams with reverse 
light quark lines are also taken into account). 

In our work, all outgoing particles are taken to be on- 
mass-shell, k\ = k\ = fc| = and P 2 = Af|. Here 
the light quark q (with momentum fe) and antiquark q 
(with momentum k^) are taken to be massless, as for the 
photon (with momentum k\). We can see explicitly that 
there exist both collinear and soft singularities in the loop 
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(c) 



FIG. 1: Feynman diagrams for T — > 7 + q + q with the color singlet bb 



integrals within the region where the light quark and an- 
tiquark interact with two gluons and are connected with 
a massless quark propagator. In general, these soft and 
collinear mass singularities may be separated by artifi- 
cially introducing a quark mass and a gluon mass. In the 
present case of dijet production, in order to define two 
jets, the jet energy has to be higher than certain energy 
scale, and the two jets have to be well separated by in- 
troducing certain value for the angle between q and q. 
With these two infrared cuts, the sensitivity of soft and 
collinear singularities will become much weaker. 

We use the NRQCD description for the bb bound state. 
In the static approximation, in which the relative motion 
of the bottom quarks in the bound state is neglected, the 
decay amplitude factorizes into a short distance ampli- 
tude M{bb — > jqq), with bb in the color-singlet state with 
zero relative momentum, and a long distance matrix ele- 
ment which is related to 71(0), the wave function of T at 
the origin. 

Within the static approximation, we write down the 
amplitude of Fig. 1(a) (for other related processes, see, 
e.g., Ref.0) 



C F g*eQ TZ{0) 



M 

8y/N c nM 2p ■ ki 



.(*) 

; A,fci 



rf 4 



<1 



UillAlvVj (3) 

m)7 M (^ — $1 + m)7 A l 



(2tt)4 



q 2 {q - h) 2 {q + h) 2 ({q + fo- pf 



where p = P/2, m = M/2, i and j are color indexes, g s 
and e are the strong and electromagnetic couplings re- 
spectively, and Q is the magnitude of the bottom quark 
charge in units of e. Charge conjugation invariance im- 
plies that the Feynman graphs are symmetric under re- 
version of the fermion flow. All six amplitudes con- 
tributing to the leading-order cross section are propor- 
tional to the same color factor l/(2^/W c ) S ab , and the 
5 a b combined with the color matrix of the right hand 



fermion flow contribute a factor Cf together with a color 
unit matrix. Gauge invariance ensures that we may 
sum over the photon spin by employing the substitu- 
tions J2 £ fj- £ t — ~9fiv To sum over the T spin, we use 
£ e p T e { 4 )a = -gP° + PPP^/M 2 . 

In our calculation, we choose two parameters to de- 
scribe the jets: one is the minimum energy cut of the jet, 
which is, say, about 1 GeV, and the other is the angle 
between the two jet axes, which is, say, about 36° (these 
values are just for illustration). We choose the following 
region: 



x > s , \z\ < zq, 



(4) 



where x is the energy fraction of the jet (k = mi), 
and z=cos8, 9 is angle between two jet axes. 

Since we have taken some steps to separate the soft and 
collinear divergences, we can complete the phase space 
integral in 4-dimension. The standard form of three body 
phase space in 4-dimension is: 



d 3 h 



»=i 1 y ' i 



With the parameters 

2Pki 



M 2 



, ^ x i 2, 



we can get 

rf$ 3 = 



M 2 



2(4tt) 



■^1^(2-^x0. 



With the following constrains: 

so < x 2 < 1, s <x 3 <1, I cos 6»| < z , 



(5) 



(0) 



(7) 



(8) 
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we can get the integration region of phase space: — zq < 
cos 9 < z for So < X2 < (2 — 2s Q + r)/(2 + s z — s ); ~z Q 
< cos6> < (2- 2s + r)/(s x 1 ) + 1 - 2/s for (2 - 2s + 
r)j (2 + s z - s ) < x 2 < (2 - 2s + r)/(2 - s ^o - s ), 
where r =m 2 / m 2 . Here for numerical estimate we choose 
so = 0.2 and zo = 0.8, corresponding to the jet energy 
cut 0.95 GeV and the angle 36° between q and q. 

In our calculation, we define (O^ 1S '( 3 Si)) ~ 
|^(0)| 2 [ll, and take m » M/2 = 4.73 GeV, and 

a s = 0.184, |7^(0)| 2 = 7.12 GeV 3 , (9) 

which are extracted from the observed decay widths of 
T — > 3g and T — > e + e _ with theoretical expressions for 
these decay widths including next to leading order QCD 
radiative corrections. We set the gluon mass as 10 -6 - 
10~ 8 GeV, together with the light quark mass m q as 
10 -2 - 10 -3 GeV, and find the numerical results are stable 
as these masses tend to zero. Finally we find the branch- 
ing ratio for the dijet production in the color-singlet sec- 
tor to be 

Br(T -> jqq) = 8.2 x 10~ 4 . (10) 

III. COLOR OCTET CONTRIBUTIONS 

Next we take the color octet contributions into con- 
sideration. With the velocity expansion in NRQCD the 



wave function of T contains high Fock states with color- 
octet bb pair and associated soft gluons . We can see 
that the color-octet bb operators Og( 3 Si), OsC'Sq), and 
Os( 3 Pj) will contribute to the process T — * -/qq. The 
corresponding diagrams are shown in Fig. 2. 

The color octet operators e.g. Og( 3 Si) and Os( 1 5'o) 
are given by 

OsCSo) = ^T a X ^T a i>. (11) 

Differing from the color singlet case, the color projection 
of the octet is 

(3i;3j| 8a) = y/2T*. (12) 
The spin projection operator e.g. Hqq{P,M) is [141 ]: 



n 00 (P,M) = —^={P + M) 75 . (13) 



We can write down the decay amplitudes of the color 
octet bb as shown in Fig. 2 (a), (b), and (c) respectively: 



M a 

M b 
M c 



iV2g 2 s e Q b Tr[T a T b ]e[*> (h 



Tr[n 15 jP,M)7' I 7 Q 7 > 



(k 2 + k 3 ) 2 [(p-k 1 ^-m 2 } 

(k 2 + k 3 ) 2 [(p- ki) 2 - m 2 
2 kf 



(p — ki) 2 — TTV 



Tr[IL ls MM)^^-hh X }], 



(14) 



where Q q is the light quark charge. Then we can express 
the decay width as functions of the jet energy cut Sq 
and the jet angle cut zq. With so=0-2 (corresponding to 
the jet energy cut 0.95 GeV) fixed, the decay width (in 
arbitrary units) as functions of zq are shown in Figs. 4- 
6. We see that the decay width increases nearly linearly 
as zq increases when zq < 0.8, whereas when zq > 0.8 
the decay width will increase rapidly. This may imply 
that taking the jet angle cut to be zq < 0.8 (0 > 36°) 
is a reasonable choice to avoid the sensitivity of infrared 
singularities. 

Choosing a s (Mr)=0.184 and a=l/137, and includ- 
ing color-singlet and various color-octet contributions, we 



can write down the branching ratio as 

Br{T^ iq q) = 2.4 x 10- 4 (T|d( 3 5i)|T) 

+ 0.061 x (T|G 8 ( 3 5i)|T) 

+ 0.084 x (T\OsCS )W 

+ 0.043 x (T|O 8 ( 3 P )|T) (15) 

For numerical results of the decay branching ratio, the 
main uncertainty comes from the estimates of the color- 
octet matrix elements. The most direct way is to use the 
velocity scaling rules. However, it is shown in Refs.flrl 
[T7] ] that by using renormalization group equations the ob- 
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(a) 




FIG. 2: Feynman diagrams with the color octet bb . 
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FIG. 3: Decay width (in arbitrary units) as a function of the 
angle cut from the color octet 3 Si operator. 
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FIG. 5: Decay width (in arbitrary units) as a function of the 
angle cut from the color octet 3 Pj operator. 



FIG. 4: Decay width (in arbitrary units) as a function of the 
angle cut from the color octet 1 So operator. 



tained color octet matrix elements are smaller. Nonethe- 
less, in our calculation we still choose velocity scaling 
rules with v 2 = 0.08 for T[l5| as a roughly estimate. An- 
other estimate can be made from the T p roduction at 
the Fermilab Tevatron using NRQCD [LJ, and the ex- 
tracted matrix elements are shown in Tabled (for more 
discussions on the color octet matrix elements, see e.g. 

Refs.fllGii). 

As the first estimate, we set the color octet 1 5o and 3 Po 
matrix elements to be zero, and use the following velocity 
scaling rule relation for the 3 Si matrix element (l5j: 




TABLE I: Color octet matrix elements for T (in units of 10 2 
GeV 3 ), taken from Table V in Ref.[H||. 

and we then get the color octet 3 Si contribution to be 

Br(T -> jqq) = 1.33 x 10~ 3 , (17) 

which is larger than the color singlet contribution by a 
factor of 1.6. 

Next, we use the results from the analysis of T pro- 
duction at the Fermilab Tevatron, and we adopt the 
assumption that the production matrix elements are 3 
times larger than the annihilation matrix elements just 
like the case of color singlet matrix elements. Note 
that the color octet matrix element (Oj^ 1S ' ) ( 1 S'o)) and 
-tt (Og 1S '( 3 Po)) can not be determined independently, 
and the values shown in Tabled from [l3| are a sort of 
simplification. We use the central values in the first col- 
umn (the 3 Po matrix element is zero), and in the second 
column (the So matrix element is zero) respectively to 
calculate the decay branching ratio, and the results are: 



(T|0 8 ( 3 Si)|T) = u 4 (T|Oi( 3 Si)|T). 



(16) 



Br(T -> jqq) = 4.21 x 10~ 3 , 
Br(T -> -yqq) = 9.31 x 10~ 3 . 



(18) 
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the result very mildly, and change of the angle cut brings 
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FIG. 6: Decay branching ratio angular distributions. The 
dotted line is for the color singlet contribution, dashed line 
for 3 Si and 3 Pj color octet contributions, and solid line for 
3 Si and 1 So color octet contributions. 

The first branching ratio is contributed from Os( 3 Si) and 
OaCSa), and the second is from C 8 ( 3 5i) and O s ( 3 Pj). 
These branching ratios are really big! To compare with 
the contribution from color singlet, we draw two group 
curves based on the latter calculation. Fig. 6 shows angu- 
lar distributions of the decay branching ratio from differ- 
ent channel contributions. We see that both in the large 
angle region and the small angle region, color octet gives 
greater contribution than singlet does. If experimentally 
the jet angular distribution can be measured then we 
might find useful hints of color octet contributions. 

Fig. 7 shows the branching ratio angular distributions 
for each individual color octet operators. Here we choose 
(0j (1S) ( 3 Si))=2.O x 1CT 2 GeV 3 . The values for other 
two operators are from Tabled. These curves are also 
useful in distinguishing between contributions of different 
color octet operators. 

IV. SUMMARY 

In this paper, we studied two mechanisms i.e. the 
color-singlet and color-octet mechanisms that contribute 
to the process T — > jqq. With the help of FeynCalc 
(20| and LoopTools[21|, we calculated numerically the 
radiative decay rates by introducing the energy cut for 
the q and q jets, and the angle cut between q and q 
jets. These cuts correspond to a minimum jet energy 
of about lGeV and the angle between two jets larger 
than 9 — 36.8°. Experimentally, these cuts are necessary 
for the two jets measurement. Theoretically, introducing 
these cuts greatly reduces the sensitivity to the soft and 
collinear singularities in the loop integrals. Furthermore, 
from our studies we find that change of energy cut affects 



FIG. 7: Angular distributions for different color octet opera- 
tors. 



a nearly linear effects in the region \cos9\ < 0.8. This 
implies that in this region those large logarithms which 
are introduced by regularizing the soft and collinear di- 
vergences only give small contributions and may then be 
ignored. Therefore, our numerical results should be rea- 
sonable. 

Based on our calculation in NRQCD, we predict the 
branching ratio of T — > jqq to be 8.2 x 1CP 4 if contribu- 
tions from the color singlet bb are considered only . 

On the other hand, the color octet bb operators may 
play an important role in this process. When we take 
them into consideration, the branching ratio could be 
enhanced by about an order of magnitude [see Eq.(18)]. 
This result is very encouraging. In this sense, the pro- 
cess T — > jqq may provide an interesting test for the 
color octet mechanism in NRQCD. However, large un- 
certainties in the estimate of color-octet matrix elements 
would prevent us from making definite predictions for 
this process. A better understanding of the color octet 
matrix elements in the future will certainly improve our 
predictions. 

We hope that the measurement on T — > jqq will clarify 
the theoretical issues presented in this paper. 
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